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The velocity Ures of resonant tunneling electrons in finite periodic structures is analytically cal- 
culated in two ways. The first method is based on the fact that a transmission of unity leads to a 
coincidence of all still competing tunneling time definitions. Thus, having an indisputable resonant 
tunneling time Tros, we apply the natural definition Wres = i/TVes to calculate the velocity. For 
the second method we combine Bloch's theorem with the transfer matrix approach to decompose 
the wave function into two Bloch waves. Then the expectation value of the velocity is calculated. 
Both different approaches lead to the same result, showing their physical equivalence. The obtained 
resonant tunneling velocity v^cs is smaller or equal to the group velocity times the magnitude of the 
complex transmission amplitude of the unit cell. Only at energies where the unit cell of the periodic 
structure has a transmission of unity «ros equals the group velocity. Numerical calculations for a 
GaAs/AlGaAs superlattice are performed. For typical parameters the resonant velocity is below 
one third of the group velocity. 
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I. INTRODUCTION 

There has been an ongoing debate about the time 
an electron spends when it passes through a classically 
forbidden region (e.g. a rectangular barrier) for many 
decades. Despite a number of review articlesi^^ and many 
papers up to now there still exist different definitions of 
tunneling time. Most of the studies have been performed 
in one of the following frameworks: (i) wave packet 
analysisi*^*^i^, (ii) dynamic paths including Feynman 
path£^-8>9, (iii) physical clocksi2iiLi2^aMiiS*iSiii, (iv) 
flux-flux correlation functionsi^, (v) theory of weak 
measurementsi^, and (vi) combinations of the former—. 

The real phase or delay time 
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where t{E) is the complex transmission coefficient as a 
function of energy, and oj is the angular frequency, arises 
from a stationary-phase argument for the transmitted 
wave packet. Many of the approaches (ii) to (vi) result 
at first in complex transmission tunneling times, given 
by one of the following expressions: 
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where d/SV{x) denotes the functional derivative with re- 
spect to the potential V{x), and r{E) and r'{E) are the 
refiection amplitudes for particles coming from the left 



and right side, respectively. The corresponding reflection 
times T^(E) are given by the substitution of t{E) by r{E) 
in the left equalities in Eqs. |(5J) to l@J. In the framework 
of physical clockg-^^d^i^^d^d'^d^i^^i^'^ the resulting times 
have been the real or imaginary part or the absolute value 
of one of the times in ^ and . 

On the other hand, it was shown that the simulta- 
neous process of (a) the determination whether a par- 
ticle is transmitted and (b) if so, how long it took to 
traverse the barrier, corresponds to two non-commuting 
observables2i. There was also some direct evidence that 
the imaginary part of the tunneling time results from 
the backaction on the particle due to the measurement 
process . 

In contrast, the time an electron spends under the bar- 
rier, either finally reflected by or transmitted through the 
barrier, is consistently given by the dwell tim^^, also 
called sojourn time, which is defined as the ratio of the 
number of particles within the barrier (extending from a 
to b) to the incident fiux^^ : 



td{E) = — I m^dx. 



(5) 



A Hcrmitian sojourn time operator exists^^, which shows 
that this time is measurable. 

Again there has been no general agreement whether a 
relation of the following form must hold: 



Tn{E) = \t{E)\\T{E) + \r{E)\\R{E). 



(6) 



Based on the argument that reflection and transmission 
are mutually exclusive events, that exhaust all possibil- 
ities in the sense of Feynman, this relation served as a 
point of focus in an early reviewi. The complex tunneling 
times T^j^ fulfill Eq. Nevertheless, the arguments for 
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Eq. (|SJ have also been critisized arguing that the used ap- 
proach goes beyond Feynman's original interpretation^. 



Tunneling times in case of transmission equals unity 

Given that |t| = 1, for a certain energy, it was shown 
that phase time and dwell time are identical, not only 
for the single barrier (at energies higher than the po- 
tential energy of the barrier)^, but also for arbitrary 
structures''^. The tunneling times t^,t^ and t|.^, 
Eqs. (jSJ to (@J, also simplify to the phase time Tphase, 
Eq. for \t\ = 1 in any arbitrary structure. Thus we 
have 

\tiE')\ = 1 ^ 

rphasc(i?') =rDiE') = T^iE') = 4{E') = t'/ [E'). 

(7) 

In accordance with the real character of the phase time 
the aforementioned problem of non-commuting observ- 
ables vanishes in the case \t{E)\ — 1, since all particles 
tunnel finally through the structure; there is neither re- 
flection nor interference between reflected and transmit- 
ted particles. 

Stimulated by comments in Ref. and recent results 
from the theory of finite periodic systema^^^S^, here we 
study the tunneling time of electrons that tunnel reso- 
nantly through finite periodic systems. In contrast to 
tunneling through single barriers, periodic systems have 
a transmission probability of unity below the barrier po- 
tential at the individual transmission resonances, which 
form allowed energy bands. Due to Eq. ||7J) we can choose 
any time definition, Eqs. to lO, to calculate the tun- 
neling time at resonance. We will use the phase delay 
time, Eg. dH), and make use of the results obtained in 
Refs. and recently. 

At first view, a similiar approach to the phase time is 
the concept of the group velocity 
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which is the velocity of the envelope of a propagating 
wave packet in a medium. Here lo is the angular fre- 
quency and k is the wave number. The function uj[k) is 
normally referred to as dispersion. The solutions of the 
Schrodinger equation for a periodic potential also yield a 
(band) dispersion relation between the Bloch wave num- 
ber q and the angular frequency uj. Using E = hxj, the 
group velocity then reads: 



(9) 



In Ref. i the following relation between the tunneling 
time and the group velocity was given (for the tight bind- 
ing limit), neglecting terms due to the matching of the 
wave functions at the ends of the system: 
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FIG. 1: Periodic potential (drawn for five periods) embed- 
ded between two infinite half-spaces. Arrows denote the wave 
function (plane waves). 



Here L is the length of the periodic system. Recent nu- 
merical calculations for finite periodic systems embedded 
in regions of constant potential showed that the equal- 
ity in H10|l . does not hold in general^l. At first glance 
it might seem paradox that the group velocity and the 
phase delay time, which are both based on a wave packet 
approach, lead to confiicting results. 

The paper is organized as follows. In section ^] we in- 
troduce the transfer matrix approach, which is used in lllll 
to calculate the resonant tunneling time and the corre- 
sponding velocity. In section Hvl Bloch's theorem is used 
together with the transfer matrix to decompose the wave 
function inside the finite periodic system (EPS) in two 
Bloch waves. It is shown that the velocity operator must 
have real expectation values at transmission resonances. 
Then the velocity expectation value is explicitly calcu- 
lated. In section^we note that both velocity approaches 
lead to the same result and derive an upper bound for 
the velocity. The calculations show that Eq. H10|l taken 
from Ref. Q gives a wrong estimation when the unit cell 
transmission amplitude is small. 

Further the special case when resonant tunneling velo- 
city and group velocity are identical is discussed. 

Our results are not restricted to a tight binding model, 
but are exact as long as only coherent transport is con- 
sidered. 

Finally, we apply the analytical results to a semicon- 
ductor superlattice and illustrate them with a compila- 
tion of graphs in section IVll 



II. FINITE PERIODIC SYSTEMS 

Our one-dimensional model system consists of an 
n-fold periodic structure, extending from a: = to 
X = L = nd {d is the length of the unit cell), embedded 
between two semi-infinite half spaces with zero potential 
(Fig. 111).^ 

Assuming a plane wave, exp(i/ca;), travelling from the 
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left towards the finite periodic system (FPS), the wave 
function is given by 



t(")exp[i/c(x- L)] 



X < 0, 
< X < L ^ nd, 
X > L, 

. ^^^^ 

where k — \/2mE /h is the electron wave vector in the 
semi-infinite half spaces, and r^"^ and i*^") are the com- 
plex reflection and transmission coefhcients of the n-fold 
periodic structure, respectively. 

We start by briefly reviewing some important proper- 
ties of one-dimensional finite periodic systems^^iS^. The 
wavefunctions at the left and right interface of a certain 
region, l and "^r, respectively, are related by the trans- 



fer matrix M through^S [ 
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where 
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spin, the time reversal invariance and the conservation of 
the probability density current lead to the structure^^ of 
the transfer matrix M 



M = 



a b 

b* a* 



(12) 



where additionally det AI = 1 holds {x* denotes the com- 
plex conjugate of x). In terms of the transmission and 
reflection coefficients t and r, the transfer matrix can be 
written aa^ 



M 



1/t r*/t* 
r/t l/t* 



(13) 



Since by construction the transfer matrix of a sequence 
of layers is the product of the transfer matrices of each 
layer, the transfer matrix of a potential consisting of n 
periods is the n-th power of the transfer matrix of one 
period: 



(14) 



For n > 2, a^") and 6^"-' can be expanded to^ 

a(") = aC/„_i(Re{a}) - [/„_2(Re{a}), (15) 

= fe[/„_i(Re{a}), (16) 

where Un{x) denote the Chebyshev Polynomials of the 
second kind. The transmission T^^"^ of any (field-free) 
n-fold periodic structure is given bjsMiSSiSSi^SiS^ 
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= [l + |6pC/,ti(Re{a})]-\ (17) 



Resonances with T'"' = 1 occur if and only if 6'"' — 
&[/„_! (Re{a}) = 0. This leads to the condition of the 
transmission resonance s ^ "A« 

r(") = 1 ^ Re{a} = cos(j7r/n), j = l,...,n-l. (18) 



The corresponding Bloch wave vectors are given by 
the condition Re{a} — cos qd: 



ITT 
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— -L,...,.^ (19) 
Inserting Eq. p8|l into H15|l we get at resonance 

a(")(q™'') = t(")(g;-) = {-ly. (20) 

III. VELOCITY FROM TUNNELING TIME 

Recently the phase time for a system with n periods, 
'''phase been calculated with the help of Eqs. l(T^ . 
and fT3|> S^: 



'phase — 



(21) 
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where T*^"^ is the transmission probability of the perio- 
dic structure given by Eq. p7l) . Here we are only in- 
terested in the phase time at resonance energies of the 
transmission, where the Bloch wave vector qj^'^ = jn/nd 
and r*^"' = 1. For energies where T^") = 1 the phase 
time is equal to the tunneling time as we discussed in the 
introduction and Eq. I|21|) reduces to the in-resonance 
phase timeS^ or resonant tunneling time: 



Im{a} dKe{a} 
1 - Re^{a} dE 



(22) 
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Clearly, Trcs is proportional to the number of periods n. 
Using the natural definition 

v,,s^L/t,^^1 (23) 
we get the following resonant tunneling velocity 



,1 -Rc^ia} / 9Re{a} 
= ^ ^ -Im{a} [-^^ 



(24) 



which does not depend on the number of periods n. For 
the sake of completeness we give the result in terms of 
the transmission amplitude i = 1/a of the unit cell: 
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Now it is interesting to compare this tunneling velocity 
Wres to the group velocity Vg. From the dispersion relation 
Re{a} — cos{qd) we obtain 
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Using Eq. H24|l . the resonant tunneling velocity and the 
group velocity are related by 



1 - Re^ia} 
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Note that the equations for the resonant tunneling time 
and the resonant velocity are only meaningful if they are 
evaluated at energies where the FPS has a transmission 
probability of unity. 

Nevertheless, by increasing the number of periods to- 
wards infinity the discrete curves tend to continuous ones, 
just like the dispersion relation does. 



IV. EXPECTATION VALUE OF VELOCITY 
OPERATOR 

In the following we calculate the expectation value of 
the velocity operator applied to the exact wave function 
inside the FPS. 



translation invariance of the Hamiltonian in a finite sys- 
tem, the wave functions are still exactly given by super- 
positions of Bloch waves, Eq. H28|l . 

If we neglect spin the time inversion symmetry of the 
Hamiltonian leads to Kramers degeneracy, i.e. -^(q) = 
i?(-q). 

Therefore, the two linear independent solutions of 
Schrodinger's equation for the one-dimensional perio- 
dic potential in a finite domain are given by and 
— \E'^?g(x). Having clarified this we use 
the following ansatz for the wave function in any one- 
dimensional n-fold periodic system inside an allowed 
band: 



(n) 
FPS 



(x) —a^J^'^Uq{x) exp{iqx) 

+ exp(— igx), < x < nd (30) 



Uq{x + d) = Uq{x) , U*{x)~U_q{x). 



The dimensionless coefficients ckg"-* and al_"^ are deter- 
mined by inital value conditions, i.e. by ^'(0) and ^''(0). 
In the following we choose the Uq{x) to be normalized 



A. Decomposition of 'i' into two Bloch waves 



dxu*{x)uq{x) — d/{2T:). 



(31) 



In 1929 Bloch calculated the eigenfunctions of the 
Schrodinger equation for crystal latticcs'^^. He modelled 
the lattice by a periodic potential which spans the whole 
space. Using group theory together with periodic (Born 
von Karman) boundary conditions he proved that the 
base solutions of the Schrodinger equation for a lattice 
potential, i.e. U(x) = U(x + Ri), where is any lattice 
vector, are of the form: 

*f(x) =uq(x)exp(zq-x). (28) 

Here q is in modern terms a reciprocal lattice vector and 



Uq(x) = Wq(x + Rj) 



(29) 



is a lattice periodic function. This fact is well known 
as Bloch's theorem. Since Bloch's theorem was origi- 
nally derived for an infinite domain, this infinite domain 
is sometimes considered as necessary. 

Therefore, we will start by briefly showing that Bloch's 
theorem gives also the exact wave functions in finite sys- 
tems. 

Mathematically spoken, for the case of the infinite 
periodic potential V'M? M considered by Bloch, the 
domain of the Schrodinger differential equation and its 
solutions is £> = R'^. If we reduce the domain of the 
periodic potential (and thus of the Schrodinger equation) 
to any finite subdomain Z? C R'^ and choose an arbi- 
trary potential in the domain M.^\D the base solutions of 
the differential equation inside D are not changed. The 
physical consequence is that despite the loss of the global 



B. The velocity operator has real eigenvalues at 
transmission resonances 



Any physical state is represented by a wave function 
which is an element of the Hilbert space of the square 
integrable functions Ti = L^(R^). The Hermiticity of the 
velocity (and the momentum) operator comes from the 
fact that the normalization (^|^) = 1 leads to ^(x) = 
as |x| oo. 

The states we consider, Eq. (|ll|l . are scattering states 
which do not belong to i^(R'^). Therefore, we can not 
expect that the velocity operator applied to any finite 
subdomain is Hermitian. Nevertheless, we show that 
the expectation values of the velocity operator inside the 
periodic structure are real for resonant tunneling states. 

The calculation is restricted to the one-dimensional re- 
gion between x = and x = L inside the FPS. The ex- 
pectation value of the velocity of any state \^) is given by 
(v) ~ (vE'l'OI^') / (vE'lvE'). In our case the numerator can be 

denoted as ('I'pp^gl-OPppsl'J'FPs)' '^here Pfps is the pro- 
jection operator onto the space region of the FPS, 



Ppps — 



Partial integration yields: 



dx I 
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At resonances with | = 1 the term 
and we obtain: 



(n) 12 
FPSi 



(^^P^SI^^FPSI^P^S) = (^FPsl^^FPSl^-FPS/ ' 

i.e., the velocity has a real value. 

C. Expectation value of the velocity at resonance 

The expectation value of the velocity in the rt-periodic 
system at resonance is given by 



I* 
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vanishes 



(34) 
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With the help of Eqs. (jSO)) and the numerator is 
given by 



'fps|v^fps|*fps/ 2tt 
The denominator of (|35|l is given by 
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(36) 
■ (37) 



A detailed derivation of both equations is given in ap- 
pendix ^ and respectively. We stress that these sim- 
ple expressions hold only in the resonant case. From Eqs. 
(|35|l to (|37|l we obtain the velocity expectation value at 



resonance: 
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An interpretation of this formula will be given in the next 
section. 



D. Determination of Bloch coefficients 



(n) / (n) 



For 



Next we have to determine the ratio 

the calculation of the expectation value it would be suf- 
ficient to consider only the wave function at resonance, 
nevertheless we will derive a more general result which is 
valid for any q value. 

We rewrite the ansatz for the wave function H3UI) in a 
reduced form 



^FPs(^) ~ Uq{x) e:>q){iqx) + a!^'J'u*(x) ex.p{~iqx), 



Uq{x) = a("^Uq(x), 



~(n) ( 



(39a) 
(39b) 

(39c) 



By formally replacing q with —q and taking the complex 
conjugate in Eq. (|39cl) 



OLq = 1 



(39d) 



is obtained. 

As stated above, the remaining coefficient a^^^ is de- 
termined by the initial values ^1/(0) and ^'(0), obtained 
from the continuity of the wave function and the proba- 
bility current density across the boundary at a; = 0. To 
simplify the algebra, we make use of the fact that these 
continuity conditions are inherently incorporated in the 
transfer matrix approach. Therefore, to determine 
we use the continuity of the wavefunction at a; = and 

in) 

two different representations of ^'pps at x = d. This way 
the derivative of the wave function is not needed in an 
explicit form.'^^ The continuity of the wavefunction 
at a; = leads to 



l + r(") ={t^(0) + 5L"j7i;(0) 



(40) 



With Eq. (|39a|l the periodicity of the Uq{x) gives the 
wavefunction a,t x — d: 



^'^^^^{d) = Uq{0) exp{iqd) + 4"^^t-g(0) expi^iqd). (41) 

On the other hand ^'ppg (d) can be obtained through the 
transfer matrix H12|l of the unit cell: 



(n) 
FPS 



(rf) = (i i)M-i(^^(i) 



(a-b). 



(42) 



Again, we note that the continuity of the probability cur- 
rent density is fully incorporated in this treatment in the 
reflection coefficient Equations (|40|l . and (|41|) to- 

gether with (|42|l form two non-linear equations in Uq(0) 

and . The solution for a^' is given by 



An) ^ Q*-6*-^ + r(")(a-5-e) 
a-6-^-hr(")*(a*-&*-^)' 



(43) 



where ^ = exp{iqd) is introduced. From Eqs. (|13|l to Hlt)|) 
we obtain 



b*Un-i(Re{a}) 



a(") aUn-i{Re{a}) ~ Un-2{Re{a})' 
Inserting into H43|) we get after some algebra 

An) _ a* -b* a(")* _a* -b* 



(44) 



(45) 



Consistent with i39dll . 5^",^ = 1/ (a^"^)* is fulfilled. 
In appendix lO Eq. 14511 is used to calculate Uq{x) and 
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^'ppg(a;). Further we prove in appendix ID li the interest- 
ing identity 



(46) 



where ji„ is the incident probabihty current. Eq. 14t)|) 
is used for the calculation of a^:^^ and aq"' in appendix 
lD3l 



As a consequence of (|45|l the absolute value of 
does not depend on n: 



Ct 



5-e 



(47) 



This shows that the ratio of the amplitude of the left and 
right going Bloch wave does not depend on the number 
of periods and is a property of the unit cell only. 

Using Re{a} = cosqd, |ap — |6p = 1, and^ = exp{iqd), 
the absolute value squared of a^j can be simplified to: 



a 



Im{a}+ i/l-Rc^{a} 



(48) 



Im{a}- ./l-Rc^{a} 



This implies that at the edges of the allowed band 
reaches unity: 

q = (0,7r/d) + 2p7r/d^ I Ma} I = 1 ^ = 1- (49) 



Taking 



(n) , (n) 



into consideration, insert- 



ing (gHl) into (IHHI gives finally 

2 



at 



"FPS/ 



1 



2 ^9 



1 - Re^{a} 

r^i — ' 

— Im{a} 



(50) 



However, this result is not true for the off-resonant case 
when the first equality does not hold. Nevertheless, the 
behaviour of (wpp^g) as a continuous function oi E ov q 
is of interest, since any rational multiple of t: / d can be 
obtained as ^rcs, by choosing proper values for j and n, 
see Eq. 

For the sake of completeness, in the resonant case, 
where q^ca is given by Eq. (HHJl, and i^") = ±1, Eq. ^ 
can be simplified to 



(51) 



V. IDENTITY OF TUNNELING TIME 
APPROACH AND VELOCITY EXPECTATION 
VALUE APPROACH 



Consequently, the tunneling time approach and the velo- 
city expectation value approach are physically equivalent 
at resonance. We use the term Urcs as the electron velo- 
city in resonant tunneling through a FPS in the follow- 
ing. From the first identity in Eq. H50() follows, that the 
electron velocity is bounded above by the group velocity 



(53) 



Here the absolute value occurs, since the group velocity 
can also take negative values for positive g, e.g. for bands 
where the energy maximum is at q = 0, while the reso- 
nant velocity is always positive. 



A. Upper bound for the resonant tunneling 
velocity 

In fact, by squaring the outer equality in f^ . an 
improved inequality compared to (|53|l can be derived. 
We use that |ap — |t|^^ > 1. Additionally, inside 
an allowed band Re^{a} < 1 holds. Therefore, < 
1 — Re^{a} < Im^{a}. Now making use of the simple 
inequality - < that holds for < a; < -y and < 2; 
we obtain 



1 - Re^{a} 



< 



1 



W{a} ~ W{a} + Rc^{a} 
Finally, this proofs 

Vres < \t\\Vg\. 



^ t 



(54) 



(55) 



In both formulas the equality holds for Re{a} — 0, 
i.e. for q = ■K/2d, which is in the middle of the band in 
q-space. This proves that the resonant tunneling velo- 
city can be much smaller than the group velocity, given 
a sufficient small transmission amplitude of the unit cell. 
Furthermore, towards the band edges where Re{a} ap- 
proaches ±1, the ratio Vrcs/vg vanishes, see Eq. (|50ll . 
Comparing with Eq. H10() , taken from Ref . we conclude 
that the matching of the wave functions at the ends of 
the system can reduce the velocity considerably. 



B. Special case: Identity between Vres and 



From the outer equality in Eq. H5U|) we can conclude 
that the equality between Vres and \vg\ holds if and only 
if \a\ = \t\-^ = 1, i.e. at energies E' with \t{E')\ = 1 of 
the unit cell of the periodic system: 

v,UE') = \vg{E')\ ^ \t{E')\ = 1. (56) 



Comparing Eqs. (|27|) and (|50|l . we get the important 
identity 



^FPS/ 



(52) 



A unit cell which is formed by a symmetric double barrier 
resonant tunneling structura^2i^ possesses the property 
|t(£")| = 1 at each transmission resonance energy E' = 
En- For this unit cell type the resonant tunneling velocity 
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Vres equals the magnitude of the group velocity Vg at all 
energies £"„. For energies different from the i?„, v^^s is 
smaller than \vg\. The only unit cell where ■^res equals Vg 
for all energies is the trivial unit cell with V{x) — 0, for 
< X < d, since it has \t{E)\ = 1 for all energies. 

Further we can derive conditions from Eq. (|5()|l for the 
identity between Vj^cs and |t)g|, in terms of the coefficients 
(") (») 

i;,es ^vg^ a^^l = 0, (57) 
i'rcs = -Vg ^ a(") = 0. (58) 

The resonant tunneling velocity equals the magnitude of 
the group velocity if and only if the wave function has 
only one Bloch component. 

Considering Eq. H5t)|l this can also be written as 

a(") =OVa^;' =0^ |t| = 1, (59) 
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FIG. 2: Schematic potential of a superlattice of the form 
(BA)". The number of periods n is chosen to be 6. 

For a SL we get the following expressions for the matrix 
elementsS^ 



yielding that the wave function inside a finite periodic 
potential consists of only one Bloch wave if and only if 
the transmission of the unit cell reaches unity. There- 
fore, we obtain the picture, that the Bloch wave moving 
into the left direction is built up by a coherent superposi- 
tion of the reflected parts of waves moving into the right 
direction. 

This solves also the puzzle we mentioned in the intro- 
duction, namely, that the phase delay time and the group 
velocity both originate from a wave packet analysis, but 
lead to different results. Often the group velocity is used 
when in good approximation no reflection occurs inside 
the medium. In contrast, here, in the case that \t\ ^ 1, 
there is reflection in each period. These reflections are 
the origin of the reduced velocity, compared to the group 
velocity. 



VI. GAAS/ALGAAS SUPERLATTICE 

All presented results are valid for arbitrary unit cells 
unless otherwise noted. In this section the results are ap- 
plied to a periodic semiconductor heterostructure. Due 
to its widespread use we choose a GaAs/AlGaAs super- 
lattice (SL). Denoting the barriers (AlGaAs) with B and 
the well regions (GaAs) with A, we can give a short nota- 
tion for the potential with n periods: (BA)". Assuming 
a stepwise constant potential function in each layer the 
transfer matrix elements can be calculated analytically. 
To keep the focus on the main topic, we will use the 
following simplifications: (i) the effective mass to* is ap- 
proximated to be the same in GaAs and in AlGaAs, and 
(ii) the effective mass does not depend on the energy. 
This simplifications can be avoided but this is left for fu- 
ture work. We present results for the lowest miniband. 
The analogous calculations for higher minibands will be 
discussed elsewhere. 



R.e{a5L} = cosh Kifc cos /cL^ — C2 sinh KLf, sin /cL^ , 

(60a) 

Imjasi} = — cosh sin /cL^ — C2 sinh Kit cos /cL^ , 

(60b) 

bsL = «ci sinh exp /cL^. (60c) 

Here Lf, and are the thicknesses of the barrier and 
the well layers, 

are the decaying electron wave vector in the barrier layer 
and the electron wave vector in the well layer, respec- 
tively and 

Ci_2 = 2^^^^^ ^■^^^)- 

The energy and the conduction band offset between GaAs 
and AlGaAs are denoted E and VJ,, respectively. 

The following SL parameters are choosen to be iden- 
tical or nearly identical to systems which were studied 
experimentally and theoretically in the past2S*2L24i^. 

For the barrier and the well width we choose Lf, = 
2.5 nm and = 6.5 nm, resp. The barriers consist of 
Al.3Ga.7As, the wells of GaAs. The number of periods 
is n = 6. A schematic of the resulting conduction band 
is drawn in Fig. [3 The calculations are performed for 
the F— valley of the conduction band for a temperature 
of 4 K. As the effective electron masses in both, GaAs 
and Al.3Ga.7As, we choose to* — 0.072too, where toq is 
the free electron mass. The conduction band offset is 
Vb = 288 meV. 

Because of the strong interdependencies between the 
dispersion relation, transmission probability, group and 
resonant tunneling velocity, and the ratio of the ampli- 
tudes of the +q and —q Bloch waves, all are shown to- 
gether in Fig.O All variables are plotted on the left as a 
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FIG. 3: Shown are relevant physical properties of superlattices. In the left column all functions are plotted vs. energy E, in the 
right column vs. Bloch wave vector q. Diamonds mark the positions of the resonant levels for n = 6 periods (band structure 
given in Fig. |5J . By changing the number of periods n, the lines connecting the diamonds are formed. The dispersion relation, 
q{E) (and E{q)), the transmission probability, T''^', the normalized coefficient of the backwards propagating Bloch wave, Iq^^^'I, 
the group velocity, Vg, and the resonant tunneling velocity Uros, and the ratio of resonant tunneling and group velocity of the 
periodic system together with the magnitude of the unit cell transmission amplitude, \t\, are shown. Details are given in the 
text. 
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function of the energy E and on the right as a function 
of the Bloch wave vector q. 

For the first miniband the resonance energies are given 
by the smallest solutions Ej of Eqs. H18|l and (|60a|l . i.e., 
by solving the transcendental equation 

cos{jn/n) — cosh kLi, cos kL^ — C2 sinh nLf, sin kLw 

(61) 

The diamonds mark the values at the resonance ener- 
gies for the system with n = 6 periods. Changing the 
number of periods n results in values that lie on the lines 
which connect the diamonds. In the limit n — > cx3 one 
would also obtain the continuous lines. 

Figure|31(a) shows q{E), and Fig. O (b) shows its in- 
verse, the dispersion relation E{q). The dispersion rela- 
tion is obtained from Rejasi} = cos qd, where Rejasi} 
is given by Eq. l|()Ua|l . 

In Fig. 131(c) and (d), the transmission probability, T^^\ 
Eq. lO, is plotted. 

From the plots of E{q) and T^^\q) the equidistant be- 
haviour of the resonant levels, Eq. H19(l . in the q-space can 
be nicely seen. Of course, due to the nonlinear dispersion 
E{q), this behaviour is lost in energy-space. 

Next, in Fig.|21(e) and (f), the magnitude of the ampli- 
tude ratio of the —q and +q Bloch waves, |<3f"j |, Eq. (|47|l . 
is shown. At the miniband edges this ratio reaches ex- 
actly unity, and it drops to a minimum around the middle 
of the miniband. 

Figure |21 (g) and (h) show the group velocity, Vg, 
Eq. H2t)|) (solid line), and the resonant tunneling velo- 
city, Uics, Eq. (|24|l (dashed line). For the n = 6 case, 
the group velocity reaches its maximum in the fourth 
resonance, while the resonant tunneling velocity reaches 
its maximum in the third resonance. This shows again 
the fundamental difference between both velocities. It 
is interesting to note that the slope of v^csiq) is zero at 
q = 0, tt/c?, while the slope of Vg has its maximum value 
there. 

In Fig. (i) and (j), the ratio v^cs/vg and the mag- 
nitude of the unit cell transmission amplitude, are 
plotted. The velocity ratio is always smaller than the 
single cell transmission, except for q ~ -k jld when both 
are equal, see Eq. H55|l . For the given parameters the 
velocity ratio is below one third. 

VII. CONCLUSIONS 

We calculated the resonant tunneling electron velocity 
in finite periodic structures embedded in regions of con- 
stant potential in two different ways and proved their 
identity. 

The first method was based on the fact that = 1 
leads to a coincidence of all tunneling time definitions, 
turning up in the literature. We used the phase time 
Tphaso — hd aTgt{E)/dE and the natural definition v = 



L/t to calculate the velocity in terms of the unit cell 
transfer matrix elements and the unit cell transmission 
amplitude, respectively, yielding Eqs. and (|2SJl. 

Bloch's theorem was combined with the transfer ma- 
trix approach to separate the wave function into two 
Bloch waves, which propagate in opposite directions, 
Eq. together with either Eq. (g3 or Eq. After 

proving that the velocity operator is real at resonance, we 
calculated its expectation value, Eq. H5U|) . as the second 
method. 

Both results are completely identical, showing the 
physical equivalence of the two approaches. The reso- 
nant tunneling velocity is in any case smaller or equal 
to the group velocity. In addition, the resonant tunnel- 
ing velocity is smaller than or equal to the product of 
the group velocity and the magnitude of the transmis- 
sion amplitude of the unit cell. Thus for unit cells with 
a small transmission amplitude both velocities can differ 
considerably. We discussed that the Bloch wave moving 
in the opposite direction of the incident electrons is due 
to reflections inside all unit cells. These reflections are 
the reason for the reduced velocity compared to the group 
velocity. At energies, where the unit cell has a transmis- 
sion of unity, only one Bloch wave remains. Consequently 
the resonant tunneling velocity equals the group velocity. 

The group velocity is often used as the speed of the 
electrons inside a superlattice. Our result points out that 
taking into account the boundary conditions reduces the 
actual velocity considerably. 

In this context the often used relation "mean free 
path equals group velocity times scattering time" seems 
questionable. We suggest that in mesoscopic structures, 
where the mean free path exceeds the length of the perio- 
dic structure it should be replaced by " mean free path 
equals resonant tunneling velocity multiplied by scatter- 
ing time^^ . The presented results support the numerical 
findings of Ref. 11 

We used electron waves throughout this paper. In ana- 
logy the main results hold for the propagation of light 
waves, of phonons waves and so on in periodic structures. 
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APPENDIX A: EVALUATION OF 

Expanding the integral with the help of Eqs. (|30|l and 
(|31|l first leads to 
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-,(«) 



(vI/f|«FFPs|vI/f) 



FPS 



^ [ 4")aL",'* f dxu,ix)e'''^^ {u,{x)e"'^) + c.c. 



(Al) 



In resonance the Bloch wave vector q is given by (|19|l . the normahzation pi() leads finally to Eq. H36|l . 
i.e., q^^^ = jir/L. For these q^°^ the last integral vanishes: 



(A2) 

due to the periodicity of both, the Uq(x) and the complex 
exponential function. Using the facti^ that 



APPENDIX B: EVALUATION OF 

Expanding the integral with the help of Eqs. and 
(ED) leads to 



Hp^^fpsI^^p^s) 



(") \ _ 



v(") 



(*",|PFPs|*f^,) 



dxul{x)e'^'^'' +C.C. 



(Bl) 



Due to the periodicity of Uq{x) and considering that 
q^^ = jir/nd, the last integral vanishes: 



dx Uq{x)e 



i2qx 



'"-1 \ 

J2exp{i2qld)] u2(x)e*29^dx 
[exp{t2TTj/n)f / u2(a;)e*2?^da; = 0, 

1=0 

(B2) 



where the last identity is due to the vanishing of the 
geometric sum. Then together with the normalization 
(EUl, we end up with Eq. lli^ . 



APPENDIX C: SIMPLIFIED CALCULATION OF 
THE WAVEFUNCTION IN THE FPS 



As a supplementary result, once a^^^ is calculated from 

(|45|l . the values of ^'pps(a;) in one period, e.g. obtained 
with the transfer matrix approach, allow to calculate the 
periodic function Uq{x). Solving (|39a|l for Uq{x) we ob- 



tain: 



Uq{x) ^ : TZm-^ cxp{-iqx). (CI) 



1 - (5 



Using (|39a|l and (|45|l . the wavefunction in the entire FPS 
can be obtained easily. 



APPENDIX D: DERIVATION OF AN IDENTITY 
AND CALCULATION OF a'g"^ AND 

1. Identity - first method 

From Eqs. ^ and i.e. 



' res 



,(") 



^[Tp^sI^-PfpsI^^p^s) 



(») 



<SI^FPS|*^;^S 



(n) 



(Dl) 



and Eqs. lO and Q), i.e. 



FPSl^FPS|*FPS 



(n) 



(D2) 



we derive the interesting identity, valid at resonance. 



1 



Jin 



(D3) 
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where the incoming probability current is given by jV 
hk/m. 



2. Identity - second method 



= exTp{ikx), and second for the wave function at 
resonance inside the FPS, i.e. ^'p'^^g(a;), both over the in- 
tervaU [0, L]. Then Eqs. 0^ and yield also identity 

inr 



This equation can also be derived directly from the 
conservation of electrical charge in the form 



(D4) 



where j is the one-dimensional probability current given 
by 



j{x,t) = Re 



ih d , 
** ( TT I * 

TO OX 



(D5) 



Since in our case and j do not depend on time, Eq. 
(|D4p simplifies to 



j{x, t) = const. 



(D6) 



3. Calculation of a,"' and a|_^ 



„(") 



Making use of Eq. (|D3|I and (|39a|l we can cal- 

culate aq"-* and a^j. Choosing the (arbitrary) phase of 
ag"\ so that a^"'' is real and positive, we obtain: 



v(") 



1/2 



^ — q ^ — q ^q ■> 



(D7) 
(D8) 



Now we integrate the constant probability current den- 
sity, Eq. (|D5|I . first for the wave function to the left, i.e. 



where is given by Eq. (|45|l . 
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